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ABSTRACT 

We measured the angular diameter of the exoplanet host star e Eridani using the 
Navy Optical Interferometer. We determined its physical radius, effective temperature, 
and mass by combining our measurement with the star's par allax, photometry from 
the literature, and the Yonsei-Yale isochrones (lYi et al.ll200ll ). respect ively. We used 



the re sulting stellar mass of 0.82ib0.05 Mq plus the mass function from iBenedict et al 
(120061 ) to calculate the planet's mass, which is 1.53ib0.22 Mjupiter- Using our new 



effective temperature, we also estimated the extent of the habitable zone for the system. 



Subject headings: optical: stars 
— techniques: interferometric 



planetary systems — stars: fundamental parameters 



1. Introduction 



Hatzes et al.l (|2000l ) discovered a planet orbiting e Eri (HD 22049, K2 V) using radial veloc- 
ity measurements that spanned 20 years. They determined the pl anet's orbit and a min imum 
mass of 0.86 Mjupiter- They assumed a stellar mass of 0.85 Mq from lDrake &: SmithI (|l993l ). who 
used spectroscopic observations to derive parameters such as effective temperature, surface gravity, 
metallicity, and luminosity and then used those parameters to calculate the star's mass. 



Di Folco et al.l (|2004l ) observed e Eri using the Very Large Telescope Interferometer (VLTI) 



during the commissioning period of the VLT Interferometer Commissioning Instrument (VINCI). 
They measured a limb-darkened angular diameter of 2.148ib0.010ib0.027 mas with a global un- 
certainty of 0.029 mas, the last of which leads to an error of 1.4%. The first two errors on the 
diameter are the statistical and systematic errors, respectively. The statistical error is based on 
the dispersion of the data while the systematic error is based on uncertainties in the calib rators' 
diame ters. Di Folco et al. then used the CESAM code version 4 evolutionary models by iMorel 
(|l997l ) to obtain various parameters for the host star, such as luminosity, effective temperature, 
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age, mass, radius, surface gravity, and metallicity. The resulting CESAM mass for e Eri was 0.83 
Mq while they list their estimated mass as 0.90ib0.10 Mq and the adopted mass as 0.84 Mq. 



Benedict et al.l (j2006l ) combined Hubble Space Telescope Fine Guidance Sensor observations 



of e Eri with ground-based astrometric and radial velocity data to produce an orbit for the system 
and the mass of the planetary companion. They assumed a mass for the central star of 0.83 Mq 
and a young age for the star, '^800 Myr. Their resulting Mpianet was 1.55±0.24 Mjupiter- 

We used the Navy Optical Interferometei0 (NOI) to measure the angular diameter of e Eri 
and then calculated its physical radius and effective temperature using the combination of our 
new measurement and other observed quantities, such as parallax, interstellar absorption, and 
bolometric correction. We provide an independent estimate of the star's mass by placing it on 
an H-R diagram using our values, which then leads to an estimate of the planet's mass. Section 
2 describes the instrument and our observing procedure. Section 3 discusses how e Eri's angular 
diameter, physical radius, and Tefj were determined, and Section 4 explores the physical implications 
of the new measurements including the mass determination and habitable zone calculation. 



2. Interferometric observations 

Observations were obtained using the NOI, an interferometer located on Anderson Mesa, AZ 



(jArmstrong et al.lll998l ). The NOI consists of two nested arrays: the four stations of the astrometric 
array (AC, AE, AW, and AN, which stand for astrometric center, east, west, and north, respectively) 
and the six stations of the imaging array, of which two stations (E6 and W7) are currently in 
operation. The baselineJl range from 16 to 79 m in length and the siderostats are 50 cm in diameter, 
of which the central 12 cm are used. The NOI observes in 16 spectral channels simultaneously over 
the 550 to 850 ^m. range. 

Each observation consisted of a 30-second coherent (on the fringe) scan in which the fringe 
contrast was measured every 2 ms, paired with an incohe rent (off the fringe) s can used to estimate 



the additive bias affecting the visibility measurements (iHummel et al.l l2003l ) . Scans were taken 
on five baselines simultaneously. Each coherent scan was averaged to 1-second data points, and 
then to a single 30-second average. The dispersion of 1-second points provided an estimate of the 
internal errors. 

e Eri was observed over three nights in 2010 using baseline lengths from 19 to 79 m (see Table 
[2]). We interleaved observations of e Eri with the calibrator star 38 Eri (HD 26574) so that every 
target was flanked by calibrator observations made as close in time as possible, which allowed us 
to convert instrumental target and calibrator visibilities to calibrated visibilities for the target. 38 



^The NOI underwent a recent name change. It was previously known as the Navy Prototype Optical Interferometer. 
^The "baseline" is the distance between the stations. 
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Eri was selected because it appears to be a single star that is significantly more unresolved on the 
baselines used than e Eri. This meant that uncertainties in the calibrator's diameter did not affect 
the target's diameter calculation as much as if the calibrator star had a substantial angular size on 
the sky. 

To establish a diameter estimate for 38 Eri and to check for excess emission from unseen 
close companions that would make it unsuitable as a calibrator, we fit a spectral energy distribu- 
tion (SEP) to pub lished UBVRIJHK ph otometry, aft e r cori verting magnitudes to fluxes using 
Colina et all (|l99d ) for UBVRI values and lCohen et al.l (|2003l l for JHK values. We used Kurucz 
model atmosphere^ with effective temperature (T^s) and surface gravity (log g) values from the 
literature to calculate the SED and determine the angular diameter that best fit the photometry. 
Figure [T] shows the SED fit. The photometry, Tefj and log g values, and the resulting diameter are 
listed in Tabled! 



3. Determination of angular diameter and Teg 



The observed quantity of an interferometer is defined as the visibility squared {V'^), which we fit 
with a model of a uniformly-illuminated disk (UD) that represents the face of the star. Diameter fits 
to V'^ were based upon the UD approximation given by = [2Ji(x) /x]^, where Ji is the first-order 
Bessel function and x = ttBOudX'^ , where B is the projected baseline at the star's position, ^ud is 
the apparent UD angul ar diameter of the star, and A is the effective wavelength of the observation 
(jShao &: Colavitalll992l ). A more realistic model of a star's disk involves limb-darkening (LD), and 
the relationship incorporating the linear LD coefficient fi\ (jHanbury Brown et al.lll974l ) is: 



1 - AtA 
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Table [2] lists the date of observation, spatial frequency (a function of wavelength and baseline), 
calibrated visibilities (V^), and errors in (aV^) for e EriQ 



The LD coefficient was obtained from Claret et al. ( 19951 ) after adopting the T^s and log g 
values required for the star. The T^s used to determine ii\ has little effect on the final ^ld: if 
Teff varies by 500 K in either direction, the resulting change in ^ld of less than 1%. Similarly, 
the final angular diameter is little affected by the choice of fix. A 10% change in the fix leads to 
a change in the measured LD diameter of less than 1%. Figure [2] shows the LD diameter (^ld) 
fit for e Eri for each night as well as the fit to all the data simultaneously. The last results in an 
angular diameter of 2.153±0.028 maJl. We then combined our measurement of the star's ^ld and 



'Available to download at "http://kurucz.cfa.liarvard.edu 



*Table [2] shows the information for a portion of one night as an example, and the full table is available on the 
electronic version of The Astrophysical Journal. 

^The fit for each individual night was within 3% of the LD diameter fit that incorporated all the data. 
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the parallax determined bv lBenedict et al.l (|2006l ) to calculate the star's physical radius: 0.74ib0.01 
Rq. Table [3] lists the measurements made here as well as other stellar parameters. 



The error for the LD diameter fit was derived using the method described in iTvcner et al 



(|2010|), who showed that a non- linear least-squares method does not sufficiently account for atmo- 
spheric effects on time scales shorter than the window between target and calibrator observations. 
They describe a bootstrap Monte Carlo method that treats the observations as groups of data points 
because the NOI collects data in scans consisting of 16 channels simultaneously. They discovered 
that when the data points were analyzed individually, a single scan's deviation from the trend had 
a large impact on the resulting diameter and error calculation. On the other hand, when they 
preserved the inherent structure of the observational data by using scans of 16 channels instead 
of individual data points, the uncertainty on the angular diameter was larger and more realistic. 
This method makes no assumptions about underlying errors due to atmospheric effects, which are 
applicable to all stars observed using ground-based instruments. Figure [3] shows the probability 
density function for e Eri's LD diameter measurement. 



Once 0LD was determined interferometrically, the Tefj was calculated using the relation 



(2) 



where -Fbol is the bolometric flux and a is the Stefan-Boltzmann constant. -Fbol was computed in 
the follow ing way: the star's V and K magnitudes were dereddened using the extinction curve de- 
scribed in lCardelli et al.l (|l989l ) and the interstellar absorption (Ay) value. The intrinsic broadband 
color {V — K) was calculated and the bolometri c correction (BC) wa s determined by interpolating 
between the [Fe/H] = 0.0, and -1.0 tables from lAlonso et al.l (119991 ). They point out that in the 
range of 6000 K > T^s > 4000 K, their BC calibration is symmetrically distributed around a ±0.10 
mag band when compared to other calibrations, so we assigned the BC an error of 0.10. -Fbol 
was determined by applying the BC and the Teg was calculated to be 5039±126 K. The star's 
luminosity (L) was also calculated using the absolute V magnitude and the BC. See Table [3] for a 
list of all these parameters. 



Because the BC is an important parameter in the T^q determination , we also derived the BC 



using 



the equation relating log(refr) and BC presented by iFlowerl (j 19961 ) ar id updated by 



Torres 



(|2010l ). We used a range of T^g values that bracketed the T^d listed for e Eri in lClem et al.l (|20o3 ) by 
450 K. The maximum difference in BC between the Flower and Alonso et al. calculations was 0.03, 
which is well within our assigned error bar of 0.10, and only changed the final T^q by a maximum 
34 K (0.7%). 



4. Results and discussion 



As a check for our result we estimated the LD diameter for e Eri using two additional methods. 
First, we created an SED fit for the star as described in Section 2 and Figure [T] shows the result. 
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The only fre e parameter for the SE P fit is the angular diameter. Second, we used the relationship 
described in iKervella et al.l (j2004l ) between the {V — K) color and log 0ld- Our measured ^ld 
is 2.153ib0.028 mas, the SED fit estimates 2.04ib0.11 mas, and the color-diameter relationship 
produces 2.04±0.82 mas. 

The main sources of errors for the three methods are uncertainties in visibilities for the inter- 
ferometric measurement, uncertainties in the comparison between observed fluxes and the model 
fluxes for a given T^s and log g for the SED estimate, and uncertainties in the K magnitude and in 
the parameters of the relation and the spread of stars around that relation for the color-diameter 
determination. All three diameters agree within their errors but our interferometric measurement 
provides an error 4 and 29 times smaller than the other methods, respectively. 

Our measurement is consistent with that derived using the FLUOR beam combi ning instrument 



on th e Center for High Angular Resolution Astronomy Array at 2.126ib0.014 mas (jPi Folco et al 



20071 ). FLUOR observes in the i^'-band at 1.94 — 2.34 /im and our measurement in visible wave- 
lengths confirms their value. This was expected because both analyses incorporate limb darkening 
in their respective wavelengths, which will compensate for any difference due to the wavelength of 
the observations. 

With our newly calculated R, Tpff, and L , we estimated the mass of the central star using 
the Yonsei-Yale isochrones (Y^ , lYi et al.l 120011 ) . Figure H] shows the tracks for several masses and 
the best fit is for a 0.82 Mq star. We combined this mass with themass function (5.9 x 10^^*^ it 
1.0 X 1O~^°M0) and inclination (30.1±3.8°) from benedict et al.l (|20od ) and calculated the planet's 
mass to be 1.53ib0.22 Mjupiter, which is consistent with that listed in Benedict et al. of 1.55±0.24 
Mjupiter- The Y^ isochrones also provided an approximate age for the system, which is on the 
order of 1 Gyr. This is a slightly older age than the results from other stud i es, which are all less 



than 1 Gyr (e.g. , ISoderblom et al. 



1991 



Henry et al. 



1996 



Song et al 



2000 



Di Folco et al. 



2004 



Saffe et al.ll2005l ). We consider our age to be consistent with those measured using other methods 
because we do not claim high accuracy in that parameter. 

We also calcu lated the size of the system's habitable zone using the following equations from 
Jones et al.l (|20od l: 



ShA^cs) = (4.190 X 10-*^ Tjj) - (2.139 x 10"^ %s) + 1.296 



and 



Sb,o{Tcs) 



(6.190 X 10"^ T"? 



(1.319 X 10"^ Tcff) +0.2341 



(3) 



(4) 



where Sb^Tf.^) and Sb^o{Tf.s) are the critical fluxes at the inner and outer boundaries in units of 
the solar constant. The inner and outer physical boundaries o in AU were then calculated using 







5*6,4 (^cff) 



and 



Sb,o{Tcs) 



(5) 
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We obtained habitabl e zone boundaries of 0.50 AU and 1.00 AU. e Eri's planet has a semimajor 



axis of 3.39ib0.36 AU (jBenedict et al.ll2006l ). so there is no chance the planet orbits anywhere near 



the habitable zone. Our values are slightly smaller than those of Jones et al. (0.572 and 1.123 AU). 
e Eri features a dust ring that was imaged at a wavelength of 850 /um u sing the Submi l limetr e 



Common User Bolometer Array at the James Clerk Maxwell Telescope by iGreaves et al.l (|l998l ) 



They found the ring has a mass between 0.01 and 0.4 MEarth; a mass comparable to th e estimated 



mass of comets orbiting in our Solar System. The disk's emission peaks at 100 //m (jPent et al 
2000) and its flux compared to the central star is on the order of 10~^, a contribution that puts it 



well outside the dynamical range of the N OI, particularly at the optical wavelengths that the NOI 



uses. At a distance of 55 AU from e Eri (ILi et al.l 120031 ) . the ring is also larger than the field of 



view of the NOI. Additionally, iDi Folco et al.l (120071 ) found no evidence of dust close to the star 



that would influence our visibility measurements. 

The Navy Optical Interferometer is a joint project of the Naval Research Laboratory and the 
U.S. Naval Observatory, in cooperation with Lowell Observatory, and is funded by the Office of 
Naval Research and the Oceanographer of the Navy. This research has made use of the SIMBAD 
database, operated at CDS, Strasbourg, France. This publication makes use of data products from 
the Two Micron All Sky Survey, which is a joint project of the University of Massachusetts and the 
Infrared Processing and Analysis Center /California Institute of Technology, funded by the National 
Aeronautics and Space Administration and the National Science Foundation. 
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Table 1. SED Input Parameters. 



Parameter 


38 Eri 


e Eri 


U magnitude 


4.53 


5.19 


B magnitude 


4.37 


4.61 


V magnitude 


4.04 


3.73 


R magnitude 


3.84 


3.19 


I magnitude 


3.68 


2.74 


J magnitude 


3.43 


2.23 


H magnitude 


3.25 


1.88 


K magnitude 


3.21 


1.78 


TeS (K) 


7079 




log g (cm s~^) 


3.66 




Av (mag) 


0.05 




Teff (K) 




5156 


log g (cm s~^) 




4.57 


E{B - V) 




0.001 


6'ld (mas) 


0.86±0.02 


2.04±0. 



Source 



Mermilliod 


(1991) 


Mermilliod 


(1991) 


Mermilliod 


^991^1 



Monet et al. (2003) 

Monet et al. (2003) 

Cutri et al. (2003) 

Cutri et al. (2003) 

Cutri et al. (2003) 

Allende Prieto fc Lambert (1999) 

Allende Prieto &: Lambert (1999) 

Neckel k Klare (1980) 

Clem et al. (2004) 

Clem et al. (2004) 

Clem et al. (2004) 

Calculated here 
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Table 2 


. Calibrated Visibilities. 










SlTia f 1 fi 1 n von 






Baseline 


Date 


JD-2450000 


(10^ cycles/radian) 






Name 


18 Nov 2010 


5518.803 


55.286 


0.399 


0.013 


AC-W7 




5518.803 


57.136 


0.380 


0.013 


AC-W7 




5518.803 


42.157 


0.679 


0.029 


AC-E6 




5518.803 


59.105 


0.408 


0.012 


AC-W7 




5518.803 


44.924 


0.643 


0.026 


AC-E6 




5518.803 


62.985 


0.338 


0.015 


AC-W7 




5518.803 


46.249 


0.616 


0.038 


AC-E6 




5518.803 


64.842 


0.341 


0.016 


AC-W7 




5518.803 


47.660 


0.607 


0.031 


AC-E6 




5518.803 


66.822 


0.268 


0.016 


AC-W7 




5518.803 


50.308 


0.510 


0.039 


AC-E6 




5518.803 


70.533 


0.274 


0.018 


AC-W7 




5518.803 


51.603 


0.525 


0.045 


AC-E6 




5518.803 


72.348 


0.241 


0.021 


AC-W7 



Note. — A portion of the table is shown here as an example. The rest of the 
table is available in the online version of The Astrophysical Journal. Approxi- 
mate baseline lengths are AC-E6 34.3 m and AC-W7 51.5 m. 
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Table 3. e Eridani Stellar Parameters. 



Parameter 


Value 


Reference 


[Fe/H] 


-0.13 ± 0.04 


Santos et al. ('2004) 


V magnitude 


3.73 ± 0.01 


Mermilliod fl991) 


K magnitude 


1.78 ± 0.29 


Cutri et al. f2003) 




0.00 


van Belle & von Braun (2009) 


BC 


0.20 ± 0.10 


Alonso et al. (1999) 


Parallax (mas) 


311.73 ± 0.11 


Benedict et al. (2006) 


Luminosity (L©) 


0.32 ± 0.03 


Calculated here 


i^BOL (10~^ erg s-i cm-2) 


98.0 ± 9.3 


Calculated here 


6'uD (mas) 


2.080 ± 0.028 (1%) 


Measured here 


6'ld (mas) 


2.153 ± 0.028 (1%) 


Measured here 


-Rlinear (-^o) 


0.74 ± 0.01 (1%) 


Calculated here 


TeS (K) 


5039 ± 126 (2%) 


Calculated here 


M (Mo) 


0.82 ± 0.05 (6%) 


Determined using model 


Age 


~1 Gyr 


Estimated using model 
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Fig. 1. — SED fits for the calibrator star 38 Eri (top panel) and e Eri (bottom panel). The diamonds 
are fluxes derived from UBVRIJHK photometry (left to right) and the solid lines are the Kurucz 
stellar models of the stars. The errors for the UBV measurements were < 1%, no errors were listed 
for RI, and the errors for JHK were 7 — 16%, which are not indicated on the plot. See Table [1] for 
the values used to create the fits. 
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Fig. 2. — e Eri LD diameter fit for each night and all the data combined. The solid lines represent 
the theoretical visibility curves for the star with the best fit ^ld, the filled circles are the calibrated 
visibilities, and the vertical lines are the measured errors. 




Fig. 3.— 
Gaussian 



Histogram representing the probability density function for e Eri's LD diameter fit. 
fit to the distribution is also shown as the solid curve. 
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Fig. 4. H R diagrams for e Eri. The solid lines arc isoclirones for the stellar masses indicated, 
the dashed line is for a mass of 0.82 Mq, and the dotted lines are for a stellar age of 1, 5, and 10 
Gyr from the bottom up. 



